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ABSTRACT 

We argue that, in general, observational studies of Be-star rotation have paid insuf- 
ficient attention to the effects of equatorial gravity darkening. We present new line- 
profile calculations that emphasize the insensitivity of line width to rotation for fast 
rotators. Coupled with a critical review of observational procedures, these calculations 
suggest that the observational parameter v sin i may systematically underestimate the 
true projected equatorial rotation velocity, v Q sini, by some tens of per cent for rapid 
rotators. A crucial implication of this work is that Be stars may be rotating much closer 
to their critical velocities than is generally supposed, bringing a range of new processes 
into contention for the elusive physical mechanism responsible for the circumstellar 
disk thought to be central to the Be phenomenon. 
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1 INTRODUCTION 

Be stars are near-main-sequence, B-type stars that have 
been observed to show Balmer-line emission. It is approach- 
ing a century and a half since the discovery of the first Be 
star iSecchilll867l) . and in recent decades compelling em- 
pirical evidence has associated the emission-line episodes 
with the formation of qua si-Kcplcrian, equatorial disks (e.g. 
iDachs et ai]ll986l : lHanuschiklll99rJ) . However, there is still 
no general agreement on the under lying physical me cha- 
nisms responsible for disk formation femith et alJl200d) . 

A phenomenologically crucial characteristic of Be stars 
is that they rotate rapidly - more rapidly than any other 
class of nondegenerate star. Virtually as soon as this was 
recognized, it was hypothesized that the equatorial rotation 
velocities, w e , may be sufficiently close to the criti cal ve- 
locity , v c , for material easily to 'leak' into a disk dStruvel 
Il93ll : here v c is the velocity at which centrifugal forces 
balance Newtonian gravity at the equator, so that the ef- 
fective gravity, g e s, vanishes). Slettebak and ot hers later 
scrut i nized this hypothesis quan titatively (e.g., ISlettebakl 
Il949l. Il96rj: ISlettebak et alJll992t . establishing the current 
canonical view that Be stars' rotation is ac tually s ignifi- 
cantly sitfccritical, wi th v c /v c ~ 0.7-0.8 (e.g., |Porterlll996t 
IChauville etafll200ll) . 

If rotation does play a direct, causal role in the Be phe- 
nomenon, then it is tempting to speculate that it does so by 
reducing the effective equatorial gravity to an extent that 
allows 'weak' processes to move material into orbit easily 
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JOwockilbOo'j) . However, if this is to be a promising avenue 
of exploration, then rotation must be much closer to crit- 
ical than is generally supposed at present: at v c /v c ~ 0.7 
the equatorial gravity is reduced by only about a factor 
2, and it requires v c /v c ~ 0.95 to reduce g e s by an order 
of magnitude. Alternatively, to launch material into orbit 
ballistically at v c /v c — 0.7 requires velocities in excess of 
100 kms -1 ; launch velocities commensurate with the sound 
speed (at which the aforementioned weak processes typically 
operate) are achieved only when a rotation rate v e /v c ~ 0.95 
is reached. Similar arguments may be brought to bear on the 
angular momentum budget needed to achieve orbit. 

In the present paper we review the bases for determin- 
ing the true projected equatorial rotation velocity, v c sini, 
and the relationship to its observational counterpart v sini 
(where i is the angle between the line of sight and the rota- 
tion axis). We argue that equatorial gravity darkening may 
plausibly have led to rotational velocities of Be stars be- 
ing systematically underestimated, and that Struve's origi- 
nal hypothesis of near-critical rotation in Be stars therefore 
has not been ruled out by existing observational studies. 



2 PREAMBLE 

For a centrally condensed star, the critical velocity of rota- 
tion is given by 

v c = y/GM./Rc = y/2GM t /3R p , (1) 

where *e' and 'p' subscripts are used throughout to denote 
equatorial and polar values, respectively (throughout this 
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paper we assume uniform angular velocity at the stellar sur- 
face). We note that, even if it were possible to measure the 
equatorial rotation velocity with complete accuracy, these 
expressions imply considerable observational uncertainty in 
establishing v c . In particular, it is customary to adopt 
masses and radii from look-up tables, on the basis of spectral 
type; and yet not only is the spectral type affected by ro- 
tation, in an aspect-dependent manner (e.g. Slettebak ct al. 
1980), but the physical relationshi p between spectral type , 
mass, and radius is also modified jMaeder fc Mevnetf 2000 ) . 
Thus systematic errors in v c /v c may easily arise from this 
source alone. 

The primary focus here, however, is on the determina- 
tion of v e sini. Almost all published surveys of rotational 
velocities in Be stars rely on one of two approaches. The 
first is to infer v e sini values from, typically, line width at 
half depth, calibratin g the results through th e standard-star 
system establishe d bvlSl ettebak et alJ <ll975l) or its extension 
to rapid rotators iSlettebakll982l) . The second is to compare 
observations with 'spun-up' versions of narrow-line stars or 
model-atmos phere flux sp ectra, typically using a convolution 
approach fcf. lGravlll992r) . We emphasize that the latter ap- 
proach is unphysical: the convolution between intrinsic line 
profile and Doppler broadening function assumes the two to 
be independent, which is clearly not the case when the stel- 
lar surface parameters (i.e., T e g & g e ff) vary with latitude 
as a result of gravity darkening. 

Both techniques imply the assumption of a single- 
valued, (approximately) linear relationship between line 
width and v c sini. As we demonstrate below, the true situ- 
ation is significantly more complicated. In particular, for a 
rigidl y rotating st ar with radiative envelope that follows the 
Ivon Zeipell (Il924l) form for gravity darkening, 

T efl oc g°J 5 , (2) 

the contribution of the cooler (and hence fainter) equato- 
rial regions to the observed spectrum may be reduced to the 
point where they are effectively invisible. A necessary con- 
sequence of such behaviour is that v sin i values determined 
without proper attention to gravity darkening will system- 
atically underestimate true projected equatorial velocities. 
This point is developed in the following sections, through 
the use of illustrative calculations and by a critical review 
of observational studies. 



3 MODELS 

In Table Q we introduce parameters for thirteen stellar mod- 
els, representative of the spectral types 2 spanning the B0-B9 
interval. The mass and polar radius of the models in Tabled 

1 Not necessarily corresponding to slowly-rotating stars 

2 We stress that here, and throughout, our use of the term 'spec- 
tral type' in reference to our models is little more than a no- 
tational convenience: although, for instance, a non-rotating star 
with the same parameters as our BO model will exhibit a spectrum 
loosely resembling a star with MK typ e BO, it should always be 
remembered that the MK scheme (e.g.. iMorgan fe Keenar]|l973t) 
is two-dimensional, and variations in a third parameter, such as 
rotation, negate any direct correspondence with physical param- 
eters. 



Table 1. Fundamental parameters for the B-type models consid- 
ered herein. 
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are taken from lPorterl il996T) . with the critical velocity eval- 
uated using equation Q. To obtain the accompanying lumi- 
nosity data we use the Warsaw-New Jersey code 3 to calcu- 
late evolutionary tracks at the masses tabulated by Porter, 
determining L„ from the epoch at which the radius coincides 
with Porter's values. 

For each of the thirteen stellar models, we construct two 
sequences, differing only in the treatment of the centrifugal 
force due to rotation. For the 'complete' set of models, the 
effects of this force are fully accounted for: the surface ge- 
ometry is that of a gravitational equipotential within the 
Roche approximat ion, and the surf ace temperature distri- 
bution follows the Ivon Zeipell lll924t) gravity-darkening law 
appropriate to radiative envelopes (cf. equation . The lu- 
minosity L* is preserved under changing model geometry by 
adjusting the polar effective temperature in accordance with 
the relation 

/ r \ 0-25 

here, a is the Stefan-Boltzmann constant, and £i is the 
surface-area w eighted gravity of the rotating star (e.g., 
ICranmerlll996ll . 

In contrast, for the 'uniform' set of models the centrifu- 
gal force is neglected, so that the star remains spherical in 
shape, and exhibits a uniform surface temperature and grav- 
ity; thes e models mimi c the traditional meth o dolog y devel- 
oped bv lCarrollI dl928Tl and lShain fc Struvd (Il929l) . which 
is the basis for the convolution approa ch described bv l Grav l 
Jl992h and used by many authors fe.g.. lHowarth et al.l l997: 
IChauville et~ai]l200ll) . 



A Details of this code are give n bvlDziembowski et all <1993t) and 
iDziembowski &: Pamvatnvkhl ll993h . the only significant differ- 
ence in the present wor k being the adoption of m ore-recent OPAL 
tabul ations for opacity Jlglesias fc Rogers 1996) and equation of 
state <Rogers et al.lll996l) . 
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Figure 1. The characteristic widths, A, of the He I A4471 A line 
profile for the B2 subtype, plotted as a function of v c sini for 
the two models and four separate inclinations considered in the 
text. The width data of the complete model (solid lines) are la- 
beled with their associated sini values; those of the uniform model 
(dotted line) all lie on the same curve. For v e sini < 0.05i> c , the 
intrinsic line broadening and blended components dominate rota- 
tional effects, and the characteristic widths should be disregarded. 



4 SPECTROSCOPIC SIMULATIONS 
4.1 Method 

For both sets of surface models, we synthesize profiles for 
He I A4471 and Mg II A4481 A (including doublet and for- 
bidden components), the lines most commonly used by ob- 
servers to determine B -star v sini value s. The calculations 
use the BRUCE code llTownsendl Il997l) . which is readily 
adapted from its original intended application in modeling 
nonradial pulsation by simply setting the pulsation ampli- 
tudes to zero. 

A large grid of precomputed conventional model spec- 
tra is required by BRUCE, to describe the angular and wave- 
length dependences of the photospheric radiative intensity 
over the range of T e ff, logg values encountered between 
pole and equator. We construct such a grid using the SYN- 
SPEC spectral-synthesis code by I. Hubeny and T. Lanz. For 
the B0-B5 models, the underlying atmosphe ric structure is 
calculated using the n on-LTE tlusty code llHubenv|ll988l : 
iHubenv fc Lanzlll995l) . under the assumption that hydro- 
gen and helium are the only sources of line opacity; for the 
later subtypes, atmospheric model s are taken fro m the line- 
blanketed LTE grids published bv lKuruczl dl993ft . 

The rotating-star spectral synthesis is undertaken over 
a range of equatorial velocities < v c < v c , and at four dif- 
fering inclinations: sini = {0.25,0.50,0.75,1.00}. We char- 
acterize the line widths of the resulting profiles by using 
the fir st zero, £i, of the Fo urier transform of the flux spec- 
trum dSmith fc Gravl Il976l) : specifically, we use A = ff 1 . 
These characteristic widths A have a nearly linear relation- 
ship with full- width at half depth, but in any case our re- 
sults are largely insensitive to the particular approach used 
to characterize line width. 

To examine the robustness of our calculations, we com- 
puted a number of variants on the basic models, including 
different choices of stellar parameters (e.g., those published 



bv lSlettebak et alJll992T) . and alternative parameterizations 
of the effects of rotation (e.g., fixing the polar temperature, 
rather than varying it in accordance with equation|HJ . These 
variants introduce only small quantitative changes in the re- 
sults. 

4.2 Results: B2 subtype 

We first consider the surface models for the B2 subtype, the 
most commo n in the sample of known Be stars (see, e.g., 
|Porteriri996l) . Fig. shows the characteristic widths of the 
He I A4471 A line as a function of v c sin i, at each of the four 
sini values considered herein, for both the uniform (A u ) and 
complete (A c ) models. We note that the four solid lines each 
span a different v c sin i range because the domain v a — [0, v c ] 
has been scaled by the four different sini values; a similar 
effect is not apparent in the dotted lines, because the A u 
curves of the uniform models are self-similar 4 , and lie atop 
one another. The latter result is a corollary of the fact that 
line widths for the uniform model follow a linear relationship 
with v B sini. 

Unlike the uniform models, the characteristic width 
curves for the gravity-darkened 'complete' models (which 
furnish a better representation of real stars) are not self- 
similar, an embodiment of departures from a linear A c - 
v e sini relationship. Such departures arise from a number of 
effects, including the variations of both line and continuum 
strength with latitude, but the dominant factor for fast rota- 
tors is the decrease in continuum emission near the equator 
resulting from von Zeipel darkening. This leads to a rela- 
tively small contribution from the fast-rotating equatorial 
regions to the spatially integrated spectrum, and in turn to 
a systematic reduction in line widths relative to the uniform 
model. Thus if v sin i values are calculated using the calibra- 
tion found for the uniform model, they may underestimate 
the true value of v a sini by ~20%, or more. Moreover, the 
near-invisibility of the equatorial regions in near-critical ro- 
tators means that increasing v e has almost no effect on line 
width for such stars. 

We illustrate this line-width redundancy in Fig. |2] for 
models with v c — 395 and 460 km s -1 (0.83 and 0.97 v c , re- 
spectively). Clearly, even with our effectively noise- free syn- 
thetic data, the helium profiles are almost indistinguishable 
in terms of their widths; there is no practicable means of 
ascertaining that the two profiles shown belong to models 
with equatorial velocities that differ by almost 20%. This 
underlines the point that the degeneracy of A c as a func- 
tion of v c must seriously compromise any attempt at devis- 
ing a v sini(A) calibration, and particularly a single- valued, 
near-linear calibration, for stars strongly affected by gravity 
darkening. 

Also shown in Fig.|2]are the Mg II A4481 A profile for the 
same v c — 395 and 460 kms -1 models. The equivalent width 
(EW) of this line increases with decreasing temperature (by 
a factor two from pole to equator in our v a — 0.95u c models), 
unlike the He I lines. This increase in relative line strength 
partly offsets the decline in continuum flux (which falls by 
factor of nearly four, from pole to equator), and as a result 

4 By this, we mean that one curve can be transformed into an- 
other by scaling both A u and v e sini by the same amount 
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Figure 2. Hel A4471 and Mgn A4481 A line profiles for equator-on 'complete' B2 models at v c = 395 kms -1 (solid lines, left-hand 
scales) and v c = 460 kms -1 (dashed lines, right-hand scales). To facilitate the comparison of line widths and shapes, the models are 
plotted on slightly different vertical scales, chosen so as to normalize the line depths. 



Table 2. Velocity deficiencies SV, as a percentage of the critical 
velocity v c , for the Hel A4471 and Mgn A448lA profiles of the 
13 model subtypes; in each case, v c /v c = 0.95 and i = 90°. 
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the Mg line widths show a somewhat greater (though still 
small) response to increasing equatorial rotation. 

4.3 Results: other subtypes 

The surface models for the other subtypes exhibit behaviour 
similar to the B2 models discussed above: with the onset of 
gravity darkening, the widths of line profiles saturate, and 
become insensitive to any further increase in the equato- 
rial velocity. To quantify succinctly the degree of this ef- 
fect, we measure the line width for complete models at 
sini = 1, Vc/vc — 0.95, and then use the uniform models 
to infer a v sini value for that width. We denote the differ- 
ence (y e sini — v sini) as a 'velocity deficiency', SV, which 
we list in Table [2] for each of the 13 subtypes considered. 

In both the Hel and Mgn lines, there is a systematic 
increase in the velocity deficiency (expressed as a fraction 
of w c ) toward later spectral types. For the most part, this 
increase arises from an enhancement of the continuum grav- 
ity darkening effect, as measured by the ratio between polar 
and equatorial continuum fluxes: in the AA4471-448lA re- 
gion appropriate to the Hel and Mgn lines, this ratio is 3.7 



for the B0 models, but grows to 7.0 for the B9 models. We 
note in passing that when bolometric fluxes are considered, 
the flux ratio is fixed at » 14.7 by dint of the fact that 
Vc/v c — 0.95, irrespective of which model is under consider- 
ation. 

The overall increase in SV/v c toward later spectral types 
is modulated by the dependence of line EW on temperature. 
As discussed in the preceding section, the EW variation of 
the Mgn line tends to counteract the continuum gravity 
darkening effect. This explains why the velocity deficiencies 
exhibited by the line are smaller than those of the Hel line, 
and also why the growth in SV/v c is seen to reverse itself 
temporarily around subtype B5. 



5 PHOTOMETRIC SIMULATIONS 
5.1 Method 

In the preceding section, we have demonstrated how grav- 
ity darkening can make a near-critical star appear, when 
observed spectroscopically, to be rotating at an ostensibly 
slower rate. In this section, we broaden our investigation by 
examining the impact of gravity darkening on the photomet- 
ric signatures of rapidly rotating stars. 

For the complete surface models (cf. Section|3J, we syn- 
thesize absolute magnitudes in the Johnson B and V bands. 
The calculations use our ROTCOLOR code, which is based 
on the BRUCE code but oriented toward photometric rather 
than spectroscopic modeling. ROTCOLOR requires an input 
grid of precomputed photometric fluxes, which we obtain 
by combi ning the four -coefficient limb-darkening data pub- 
lished bv lClarell (l200fj|) with normal-emergent absolute flux 
data made available to us by Cl aret (pe r sonal communica- 
tion); both data are based on the lKuruczl il993Tl atmosphere 
grids. To convert synthetic fluxes into absolute magnitudes, 
we use the observed va lues for the Sun, Mb = 5.48 and 
M v = 4.83 jAllenlll976l) . 

Our calculations are undertaken for five rotation veloci- 
ties v a /v c = {0.00, 0.24, 0.48, 0.71, 0.95}, and at three differ- 
ing inclinations i = {0°, 45°, 90°}. For comparison purposes, 
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Figure 3. Selected complete surface models (cf. Section plot- 
ted as points in the Johnson B — V colour-magnitude diagram, 
at five rotation velocities v c /v c = {0.00, 0.24, 0.48, 0.71, 0.95} and 
three inclinations i = {0°, 45°, 90°}. Solid lines join together the 
points from models that differ solely in their rotation velocity; 
dotted lines indicate ZAMS— TAMS evolutionary tracks that pass 
through the non-rotating models, the latter being labeled at the 
bottom of the plot with their nominal spectral types. 

we also synthesize colours for the evolutionary tracks calcu- 
lated in Section |3] each track extends from zero-age main 
sequence (ZAMS) through to terminal-age main sequence 
(TAMS, defined by the cessation of core hydrogen burning). 



5.2 Results 

In Fig. [3]we plot the results of our photometric simulations 
in a colour-magnitude diagram. The general effect of rota- 
tion is to displace models away from their zero-rotation loci, 
moving them along trajectories in the My j(B — V) plane 
that are, to a first approximation, straight lines. The dif- 
ferential rate of displacement varies in step with v c , being 
largest close to critical rotation, and so small for the slowly 
rotating models (v c — 0.24 v c ) that they can barely be dis- 
tinguished from their non-rotating neighbours. 

The most striking aspect of Fig. [3] is the strong incli- 
nation (i) dependence of the trajectories' orientation in the 
My /(B — V) plane. For the i = 0° models, rotation increases 
the observed brightness, due to beaming of radiation from 
the hot pole; however, this increase happens without any 
appreciable variation in colour, leading to trajectories that 
are close to vertical. The converse occurs when i = 90°: 
the centrifugal distention of equatorial regions produces a 
drop in the surface-averaged effective temperature, redden- 
ing the colours at approximately constant brightness, and so 
resulting in nearly horizontal trajectories. For intermediate 
i = 45° combination of brightening and reddening 

produces diagonal trajectories, coincidentally nearly paral- 
lel to those produced by the process of stellar evolution. 

The increase in brightness for the pole-on models is ap- 
proximately independent of subtype, being AMv w — (f?65 
at v e = 0.95 v c , while the reddening of the equator-on mod- 
els, at the same rotation rate, shows a modest increase from 
A(B - V) = m 032 at B0 to A(B - V) = ff?066 at B9. 
By chance, the ratios between these brightening and red- 



dening values correspond closely to the colour-magnitude 
slope of the main sequence throughout the B spectral type. 
Therefore, for an ensemble of randomly oriented B stars, the 
effect of gravity darkening due to rotation at fixed v c /v c is 
to shift the main sequence up and to the right in the colour- 
magnitude diagram, but to leave it otherwise undistorted. 
We will return to this point in Section 1531 



6 DISCUSSION 

The calculations in Section [I] are intended to underline our 
central point that line widths are a rather poor indicator 
of v c sin i at near-critical rotation. Although this basic idea 
has been previously noted in the literature iStoecklevlll968l : 
lHardorp fc Strittmatterl Il968al: ICollins fc Truaxl Il995l)~ hi 
pra ctice it has be en almost entirely ignored by observers 
(cf. lHowarthl2002D . For exam ple, in the course of writing up 
this work we discovered that IStoecklevI il96ct his Fig. 4) 
gives a diagram that is very similar to our Fig.Q while also 
addressing the implications for Be stars. However, the low 
citation rate for this paper (at the time of writing, a sin- 
gle citation 5 in the ADS database since 1975) suggests that 
the significance of this result for Be-star physics has been 
largely overlooked. In this section, therefore, we critically 
review the historical basis for the canonical view that v c /v c 
does not exceed ~0.8 for Be stars. 

6.1 Line profiles 

The basic method for determinations of v sini is, of course, 
line-profile measurements. Slettebak's (1982) observational 
study of rotation in a sample of 163 Be stars has been 
particularly influential; many authors have relied directly 
on his v sini dat a (e.g.. lPorterlll996T). or have used them 
to ca l ibrate (e.g.. lBalonalll975l: lHalbedellll99r3: ISteele et all 
[ 19991: lYudinl l200ll; lAbt et alJ 12002ft or t o validate (e.g., 
iHowarth et alJll997l: IChauville et al.ll200lft their own mea- 
surements^^ 

ISlettebakl (I1982TI used visual inspection of photographic 
spectra to estimate the line wid ths in his targets ; calib rating 
his results with standards from lSlettebak et alJ lll975l) . This 
approach assumes a single-valued correspondence between 
line width and v sini, and neglects the redundancy between 
these parameters that we have emphasized. Slettebak is uni- 
versally acknowledged as an exceptionally skilled observer, 
and we believe that his careful and conservative approach 
would have led him to adopt the minimum v sin i value con- 
sistent with his observations (rather than some arbitrarily 
larger value), which, as we have shown, will generally un- 
derestimate v c sini. 

Furthermore, the v sini values of the underpinning 
ISlettebak et alJ <ll975l) standard stars were themselve s cali- 
brated against theoretical profiles from ICollins! Jl974l) . The 

5 ICollins fc Truaxl Jl995f) cite Stoeckley's work, and also give a 
similar diagram to our Fig- □ (their Fig. 4); but we disagree with 
their conclusion that an approximate treatment of limb darkening 
is the major factor in departures from the results of the uniform 
model. Note that even our simple, uniform models incorporate 
a detailed treatment of non-linear, wavelength-dependent limb 
darkening. 
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latter's pioneering calculations required computations that 
were, by the standards of the day, sophisticated and de- 
manding . Unfortunately, subse quent results, including his 
own (cf. ICollins fc T ruax 1995) and those we present here, 
shed s ome doubt on this fundamental calibration. For ex- 
ample, ICollinsI il974T) found line broadening to be relatively 
unaffected by gravity darkening. Indeed, his paper explic- 
itly notes a "remarkably good linear correlation between 
the theoretical value for v c sini and the half- width of the 
line" - a correlation from which our gravity-darkened, com- 
plete models depart, quite significantly, in the crucial limit 
of large v c sin i (Fig. . 

Collins' remark was illustrated by calculations of the 
He I A4471 A profile for a BO- type stellar model (his Fig. 5). 
That figure indicates line half-widths at half depth that 
are systematically broader for critical rotators than for sub- 
critical rotators at the same v e sini. We find it difficult to 
devise any physical explanation for such behaviour, and sus- 
pect that it may indicate some proc edural or coding error in 
his study. Indeed, later work by, e.g.. lCollins fc Truaxl (Il995l) 
gives findings broadly consistent with our own, and also con- 
tradicts the linear correlation found by [Collins (1974J). In 
any event, our Table [5] indicates that the BO subtype ex- 
hibits the smallest velocity deficiencies of any of the stellar 
models that we have considered. Therefore, notwithstand- 
ing any putative problems with Collins' early calculations, 
we expect the data plotted in his Fig. 5 to provide the poor- 
est illustration of the degree to which gravity darkening can 
influence line broadening. 

From the point of view of both observational procedure 
and underlying calibration, we are therefore led to the con- 
jecture that t he system of Be - star v eloci ties underpinned 
by the work of lSlettebak et al.1 (Il975l) and ISlettebakl (I1982D 
may systematically underestimate v c sini values, particu- 
larly at the upper end of the rotation scale. This conjecture is 
given credence by the recent investigation bv lChauville et all 
(2001), who obtained v sini values for 116 Be stars by least- 
squares fits of synthetic spectra to observed He I A4471 A line 
profiles. As for our uniform models, Chauville et al.'s syn- 
thetic profiles were calculated under the assumption of no 
rotational distortion or gravity darkening, and they there- 
fore certainly underestimate v E sin i. However, they obtained 
v shu values in excellent overall agreement with those of 
ISlettebakl il982T) . encouraging the view that there is also 
a systematic error in the latter's results (and hence in all 
subsequent work built on them). 

6.2 Statistics 

IChauville et alJ (1200111 measured a mean projected equa- 
torial velocity ( v sini/v c ) = 0.65 for the i r sam ple of Be 
stars. Using the IChandrasekhar fc Muncbl (ll950T) relations 
between projected and intrinsic rotation velocities, they in- 
ferred from their data a mean equatorial velocity (v c /v c ) — 
0.83 ± 0.03. This extremely small apparent intrinsic scatter 
led them to conclude that "all studied Be stars rotate at 
nearly the same ratio [of] v e /v c " . 

We know of no physical reason why all Be stars should 
rotate at the same v c /v c ~ 0.8, but examination of our Fig.0 
shows that at and above v c /v c ~ 0.8 there is effectively no 
change in line width with increasing v a . It is therefore rea- 
sonable to suppose that the observations may be consistent 



with any value (or a range of values) > 0.8, and particularly 
with near-critical rotation - an obvious limiting value. 

To investigate this possibility in a heuristic way, we con- 
sider a population of complete B2-star models with equato- 
rial velocities distributed uniformly between v a = 0.93 v c 
and v e = 0.97 u c , and with rotation axes oriented ran- 
domly. Analyzing this population with He I A4471 A line 
profiles derived from our uniform models (thereby effec- 
tively mimicking Chauville et al.'s procedure), we obtain a 
mean proj ected equatorial velocity (v smil v c ) = 0.635, from 
which the lChandrasekhar fc Miinchl dl950F) relations suggest 
a mean equatorial velocity (v c /v c ) — 0.809, with a variance 
a 2 (v e /v c ) = —9 x 10 -3 . These inferred values 6 are remark- 
ably close to those obtained by Chauville et al., even though 
our synthetic population rotates at significantly more-rapid 
rates than the (v c /v c ) — 0.83 advanced by those authors. 
Of course, this simulation is not particularly realistic, but 
it nonetheless clearly establishes that statistical studies that 
do not adequately treat the effects of gravity darkening must 
inevitably lead to systematic underestimates of v e sini. 

6.3 Photometry 

The calculations in Section |^| indicate that near-critical 
rotation (by which we mean v c /v c ~ 0.95) displaces 
stars toward brighter magnitudes or redder colours, de- 
pending o n the inclination. This phenome non, first dis- 
covered by iRoxburgh fc Strittmatterl dl965l) . has been in- 
vestigated and confirmed by a number of authors (e.g., 
CollinsllT96l iFaulkner et alJll968t lllardorp fc Strittmatterl 
1968bl : iMaeder fc Pevtremannlll97fJ) . Our qualitative find- 
ings are very similar to these studies, although we obtain 
smaller changes in brightness (AMy ~ — 0^6) and colour 
(A(_B — V) ~ 0™06) than found historically (for reasons we 
discuss below). 

Evidently, if Be stars are indeed near-critical rotators, 
they should naturally occupy a brighter/redder position in 
the colour-magnitude diagram than normal B-type stars. A 
search through the relevant literature r eveals substantive ev - 
idence for such an anomalous position: IZorec fc Briotl (|l99lT> 
find their sample of Be stars to be over-luminous with re- 
spect to main sequence stars, and of the fourteen previous 
studies that they cite, eleven disclose similar behaviour. 

In the era when the observational consequences of grav- 
ity darkening were first being investigated, it was natural 
to at tribute this obs erved over-luminosity to rapid rotation 
fe.g.. lStoecklevll968h . However, many of the studies from the 
1960's period (including the majority of tho se we cite above) 
were based on interior c alculations by IRoxbu rgh et all 
lll965h . which were found bv lSanderson et alTi^Tol) to con- 
tain a serious error. Correction of this error resulted in a 
reduction of the brightness/reddening change caused by ro- 
tation (towards values commensurate with those we found in 
Section |SJ, to such an extent that rotation alone was hard- 
pushed to explain the luminosity anomaly. 

Ever since the discovery of the error bv lRoxburgh et al] 

6 Note that the variance for our synthetic population is less 
th an zero due to the marg i nal br eakdown of the assumption 
by IChandrasekhar & Munch (1950) that the v sini distribution 
arises from the spatial projection of an intrinsic distribution. 
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(1965), there has been a shift towards interpreting the 
anomaly as a consequence of the circumstella r emission (CE ) 
produced by the disks of Be stars fe.g.. lFabregat et, alll99fi|l . 
This hypothesis certainly seems promising, but it has never 
been est ablished whether CE can explain the observations in 
toto; as lZorec fc Briotl lll997f) point out, it may be the case 
that multiple mechanisms (including rapid rotation) are re- 
quired to explain the anomalous position of Be stars in the 
colour-magnitude diagram. Indeed, it is still not clear what 
th e effects are of t he intrinsic 'rotation reddening' discussed 
bv iMaedeH lll975h . which arise from the distention of the 
envelope around the equatorial regions of a rapidly rotating 
star (such that the local atmospheric structure departs from 
that predicted by plane-parallel models) . 

6.4 Future Approaches 

Having presented evidence that current measurements of Be 
star rotation are underestimates, let us briefly consider what 
approaches might be deployed in the future to address this 
problem. To reiterate, the essence of the problem is this: 
when a star rotates at a significant fraction of its critical 
velocity v c , the insensitivity of line widths toward any fur- 
ther increase in the rotation rate means that these widths 
are poor indicators of the star's true projected equatorial 
velocity v e sini. How might this degeneracy be lifted? 

One possible approach lies in leveraging the combined 
temperature/wavelength dependence of the continuu m flux. 
From far ultraviolet (FUV) observations of Be stars, iHearj 
(1976) discovered the line broadening at these wavelengths 
to be significa ntly less than a t visible wavelengths. As 
pointed out bv iHutchingsl dl976l) . this effect can be under- 
stood as a result of the different surface brightness distri- 
butions at visible and FUV wavelengths , caused by gravity 
darkening. iHutchings fc Stoecklevl (Il977li subsequently sug- 
gested that the effect might be used to obtain reliable values 
for v a /v c and i. Such an approach certainly looks promising, 
but suffers from the drawback of requiring simultaneous vis- 
ible, FUV and (if possible) infrared observations, which are 
often difficult to secure. 

An alternative technique for lifting the degeneracy lies 
in exploiting the differing temperature sensitivities of dif- 
ferent spectral lines. As we demonstrated in Sections 14. 2h 
14.31 the inverse temperature/EW dependence exhibited by 
the Mg II A4481 A line partly counteracts the effect of grav- 
ity darkening in the continuum, and lessens the impact of 
the degeneracy. In the case of an individual line such as 
this, the effect will probably be too small to be of practi- 
cal use 7 , but the effect could be exploited by simultaneous 
modelling of multiple lines. In principle, the different tem- 
perature/EW dependence of each line means that a good fit 
for all lines will occur only at a single point in parameter 
space (that is, a single combination of v e , i, M* etc.). In real- 
ity, of course, noise and other factors will limit the precision 
of this approach; nevertheless, the approach certainly ap- 
pears promising, and has alrea dy been applied with suc cess 
to rapidly-rotating O stars fsee lHowarth fc Smith 201^). 

7 We also note on more general grounds that the Mg II A4481 A 
line is a poor choice for v e sini determination, due to blending 
with the adjacent He I A4471 A line. 



The most direct method for examining rotational dis- 
tortion is through direct imaging. Once the aspect ratio 
x = R c /R p of the centrifugally-distorted star is known, its 
equatorial velocity can be calculated in the Roche approxi- 
mation from the expression 

V ~ = Wfcil, (4) 

V c V X s 

where w e note that x = 1.5 correspond s to v c /v c = 1. 
Recently. iDomiciano de Souza et alJ ll200ot) applied this ap- 
proach to the Be star a Eri; their VLTI observations in- 
dicate an aspect ratio entirely consistent with near-critical 
rotation. Clearly, the power of this approach lies in its abil- 
ity to avoid the need for complex modelling of stellar atmo- 
spheres or interiors. Unfortunately, its applicability is lim- 
ited to those stars which are resolvable via interferometry. 



7 CONCLUSION 

We have argued that the available observations, and their 
interpretation, leave open the possibility that Be stars rotate 
at or near to their critical velocities, contrary to the previous 
consensus. In questioning the canonical view that Be-star ro- 
tation is substantially sub-critical, we have been motivated 
by the issue of the fundamental physical processes respon- 
sible for the Be phenomenon. A range of phenomenological 
models have been developed to account for various aspects 
of B e-star behaviou r, some quite detailed and sophisticated 
fe.g.. lOkazakill99ll) . but the mechanisms underlying the for- 
mation of circumstellar disks remain elusive. This is largely 
due to the substantial energetic requirements for levitating 
material in a strong gravitational field. 

Near-critical rotation alleviates this problem to the 
point where known processes, such as pulsation or gas pres- 
sure, may readily provide sufficient energy and angular mo- 
mentum. Our critique argues that such near-critical rotation 
has not been ruled out by existing studies, be they spectro- 
scopic or photometric. Indeed, these studies furnish neces- 
sary evidence that the Be stars are rotating close to critical; 
we hope that this paper will provide a motivation towards 
obtaining, perhaps using one of the approaches we suggest 
in Section 16.41 sufficient evidence in support of the same 
conclusion. 
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